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ABSTRACT 

Merger shocks induce turbulence in the intra-cluster medium (ICM), and, under some 
circumstances, accelerate electrons to relativistic velocities to form so-called radio 
relics. Relics are mostly found at the periphery of galaxy clusters and appear to have 
magnetic fields at the /iG level. Here we investigate the possible origins of these mag- 
netic fields. Turbulence produced by the shock itself cannot explain the magnitude 
of these fields. However, we argue that if the turbulent pressure support in the ICM 
upstream of the merger shock is of the order of 10 to 30 percent of the total pressure 
on scales of a few times 100 kpc, then vorticity generated by compressive and baro- 
clinic effects across the shock discontinuity can lead to a sufficient amplification of the 
magnetic field. Compressional amplification can explain the large polarisation of the 
radio emission more easily than dynamo turbulent amplification. Finally, clumping of 
the ICM is shown to have a negligible effect on magnetic field amplification. 

Key words: magnetic fields - turbulence - galaxies: clusters: general - shock waves 
- methods: analytical 



1 INTRODUCTION 

In the classification of diffuse synchrotron radio emission 
in galaxy clusters, radio relics are polarised (generally at 
a level of 10 to 30%) sources with a steep spectral index 
Q > 1. They have elongated shapes, with sizes between 10^ 
and 10^ kpc, and ar e usually located at the cluster periphery 
l|Ferrari et al.ll2008l. for a review). According to the scheme 
proposed by Kempner et al.l l|2004l ). there are two classes of 
relics, the "radio ghost" cau sed by shock-induced compres- 
sion of fossil radio plasma (|Enfilin fc Gopal-Krishnal l200ll : 
lEnfilin fc Briig gcn 200j) and the "radio gisch t" caused by 
shock accelerated cosmic ray (CR) electrons (jEnfjlin et al.l 
ll998l : lHoeft fc BriiggenllioOTT l. In this paper, we will only be 
concerned with gischt relics. 

Relics are thought to be caused by shocks, triggered by 
a major merger and propagating throug h the in tr a-cluster 
medi um (IC M). In a few cases (A5 21, iGiacintucci et al 



20081: A3667, iFinoguenov et al.ll2010l : A754. iMacario et al 



20111 ) the radio emission in relics is co-located with a shock 



front, detected by X-ray observations. 

Some observations suggest that relics host relatively 
large magne tic fields, with typical values of a few /xG 
(for example. iBonafede et alJ[2009l : IFinoguenov et ahllioiol l. 
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These fields are derived using several methods, mostly 
equipartition argum ents, comparison o f hard X-ray and 
radi o emission (cf. Ferrari et all l2008l ). the spectral ag- 
ing (|Markevitch et all 20051 ) and the width of relics seen 
edge-on. The last method has been employed for a 
reUc recently observed in the cluster CIZA J2242. 8+5301 
(|van Weeren et al]l201Gl ). This relic is extreme in terms of 
its very regular morphology, the high degree of polarisation 
and the large magnetic field, between 5 and 7 /iG. 

The theoretical study of relics in cosmological simula- 
tions requires to connect the simulation output with a post- 
processing step for the computation of the radio emission. 
This exercise has been performed several times in the liter- 
ature, using different hydrodynamical codes (m esh-based or 
SPH) and assumptions for the radio emission (i Hoeft et al.l 
[2008; Battaglia ct al. 2009; Skillm an et al. 2011, ). Other re- 
lated works focus on the accelera tion mechanism of the 
CR electrons (a recent example is iKang fc Rvul 2011^■ on 
predictions for fu ture radio observations I Nuza et al.l 20121 : 
IVazza et al. 2012) and on putting constraints on the merger 
geometry (van W ccrcn ct al. 2011a). The cited studies make 
assumptions for the magnitude of the magnetic field B 
in the ICM, either bas ed on a scaling of B with density 
l|Hoeft et al] [20 08: Skill^meEiD [IqiJ), or with thermal 
energy teatta glia et al.ll2009l ). because those simulations do 
not include any MHD treatment. However, there are a num- 
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ber of observations and simulations of magnetic fields in 
clusters, for the calibration of these scalings. Here we briefiy 
review the results most relevant to our analysis. 

Observations of relics show (cf . Ivan Weeren et alll2009l ') 
that the most extended, arc-like objects are found at large 
projected distances from the cluster centre. For example, 
CIZA J2242. 8+5301 is a massiv e cluster (Lx = 6.8 x 
10** erg s"\ thence T ~ 9 keV, Ivan Weeren et all l20ldl ). 
so its virial radius would be about 3 Mpc, meaning that its 
relic is located approximately at 0.5 -Rvir. Another promi- 
nent object, the northwest relic in A3667, has a projected 
distance of about 1 -Rvir. These two cases set the distance 
range that this study investigates. The typical magnetic field 
near t he lower edge of this range, derived by IClarke et al.l 
(2001) from the statistics of Faraday rotation measure (RM) 
measurements, is 5 /iG within a distance which, given the 
X-ray luminosities of th e targets in their cl uster sample, is 
approximately 0.3 -Rvir. iGovoni et al I (|2006l ) study the RM 
in A2255 and, with their best model, derive a radial profile 
for B with values between 0.6 and 0.4 /^G in the distance 
range between 0. 5 and 1 -Rvir. Similarly , in their RM obser- 
vations of Coma, iBonafede et al.l l|2010l ) derive a profile for 
B obeying B{r) = Bt)[n(r) / noY' , where n is the gas num- 
ber density, and -Bo and no are central values. Their best 
model is for -Bo = 4.7 /iG and = 0.5, resulting in a value 
of 0.9 /iG at 0.5 -Rvir. Wit h a different technique, based on 
de-polarisation arguments. iBonafede et al.l (|2011bl) find m a 
sample of clusters on average -B ~ 2.6 /-iG within approxi- 
mately 0.3 -Rvir. 

In summary, observations support that the ambient 
magnetic field strength in the distance range relevant for 
radio relics is, within a factor of a few, 1 /iG at 0.5 -Rvir. 

These values are confirmed by cluster MHD simulations. 



tional phvsics (|Dolae et al. 


20021: Dubois & Tevssieil 20081: 


Dolag & Stasvszvnl l2009l: 


Donncrt et al.' 


20091: Xu et all 


2010: Bonafede et al.l 2011al: Xu et al. 2011 


). The theoreti- 



cal results show a marked dependence on the magnetic seed, 
either primordial or produced by an AGN outflow, but un- 
der the constraint of a central magnetic fleld of a few /iG 
they are generally in agreement with the observations. 

In this work we want to address the following questions: 
How are the apparently large magnetic fields in relics pro- 
duced? Can we explain the observed polarisation in relics? 
Our fiducial example will be the relic in the galaxy cluster 
CIZA J2242. 8-1-5301? 

An important process for magnetic field amplification, 
often invoked in the magnetogenesis of galaxy clusters , 
is the small-scale, turbulent dvnamo |Dolag et al.l l2002l : 
iBriiggen et al.ll2005l : ISubramanian et aLlbOod ). However, we 
will show that the amplification of the magnetic field, B, is 
not driven by turbulence injected by the shock. Instead, we 
argue that ICM turbulence upstream of the shock produced 
a magnetic field, that is amplified by the shock. Thus the 
upstream magnetic field can be amplified by a factor close to 
the shock compression ratio or somewhat larger if baroclinic 
effects are important. 

The role of shock-driven vorticity generation will be 
studied in Section [21 where estimates for the downstream 
magnetic field and the polarisation level will be derived. The 
efficiency of turbulent amplification and the effect of clump- 



ing will be quantified in Section [3] and [4l respectively. The 
results will be discussed and summarised in Section (5] 



2 VORTICITY JUMP ACROSS A 
PROPAGATING SHOCK 

2.1 The vorticity equation and its terms 

Hydrodynamical shocks can amplify vorticity oj which in 
the frozen-in case follo ws the s ame evolution equation as the 
m agnetic field, B (e.g. [Pavies fc W idrow 20001. According 
to lKevlahan fc Pudrit j (|2009l ') (their equation 5), the vortic- 
ity jump 5u3 across a propagating shock, along the binormal 
direction b (tangential to the shock surface) is: 



Suj ■ b = 



dMs 



l + fi dS 



+ 



1 



1 + A* 



-Mi 



d\Ml 

dS 



+ U) X U ■ S \ + ^U3 ■ b 



(1) 



where s is tangential to the shock and n is the shock-normal 
direction, so that the binormal direction is b = s x n. We 
note that the normal component of the vorticity is continu- 
ous across a shock. The normalised density jump across the 
shock is /i = p2 /pi — 1, where p2 and pi are the pre-shock and 
post-shock gas densities, respectively, Mt is the turbulence 
Mach number of the upstream flow, -Ms the Mach number of 
the shock and d/dS is the spatial derivative in the direction 
tangential to the shock surface. 

The first term on the right-hand side of equation (O, 
the so-called curvature term, is relevant at locations where 
the shock speed varies along its surface, or where this surface 
bends. The second term governs the baroclinic generation of 
vo rticity, as i t clear ly appears in the derivation performed 
bv lKevlahanI (|l997l '). but it is cast in a more general form 
than Vp X VP (e.g.. lRvu et al.|[200l '): in particular, the ex- 
pression in equation ((TJ is suited to deal with the shock 
propagation in a non-uniform flow. The third term is asso- 
ciated with flow compression at the shock, and arises from 
the conservation of angular momentum. We stress that this 
term in equation ((T)) is non-vanishing only if the upstream 
fluid is not irrotational (i.e., if it h a s a no n- vanishing vortic- 
ity); on the other hand. iKevlahanl (|l997l ) shows that even a 
straight shock can generate vorticity in an irrotational fluid 
by baroclinic effects. 

If the upstream flow is turbulent and the incoming shock 
is curved, all terms can potentially contribute to the amplifi- 
cation. The compressional contribution has the same nature 
as the magnetic amplificati on at shocks follow ing from flux 
conservation, described bv lEnfili"n" et all (|l998l '). Consider a 
magnetic flux tube with strength B\, oriented at an angle 
9\ with respect to the shock normal: this tube will be bent 
by the interaction with the shock and the field amplified, 
according to 



and 



B2 



tan62 = -Rtan^i 



-Bi (cos^ 6\ + -R^ sin^ i 



,1/2 



(2) 



(3) 



where the subscripts 1 and 2 refer to the pre-shock and 
the post-shock region, respectively, and -R = pij p\ is the 
shock compression ratio. From equation ((3} it is easy to see 
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that the magnetic field amplification is negligible when Bi 
is oriented parallel to the shock normal {9i ~ 0, thereby 
B2 — Bi), whereas B2 — RB\ when B\ is nearly paral- 
lel to the shock surface. The shock compression amplifies 
therefore preferentially the field along the shock surface, al- 
though only with a moderate efficiency namely, at most, 
-R = (7 -1-1)/ (7 — 1) for strong shocks. For a typical merger 
shock with Mach number Ms ~ 3 the Rankine-Hugoniot 
shock jump relations predict i? ~ 3, assuming a polytropic 
index 7 = 5/3. Hence, shock compression results in an am- 
plification of _B by a factor of a few. If the shock is further 
modified by cosmic rays, the compression factor can become 
even higher (cf Section [5]). 

How strong is the v orticity produced by the curvature of 
the shock? According to lTruesdelj l| 19521 ). the jump induced 
by shock curvature can be rewritten as 



compressional term 
baroclinic term, M, = 0.3 



stiock Macti number 



5uj ■ 



1 



-UsK 



(4) 



where Us is the fluid velocity tangential to the shock, in the 
shock reference frame, and K is the shock curvature. For 
a simple estimate in the merger shock case, let K ~ I/O, 
where D is the distance of the shock from the cluster cen- 
tre, set to 1.5 Mpc for similarity with CIZA J2242. 8-1-5301. 
We also approximate Us with the turbulent velocity in the 
upstream medium. To be more quantitative, the ratio of 
turbulent to thermal pressure at the length scale I can be 
written as 



-Pturb 



(0 



.(0/3 



kT/ifinip) 



(5) 



where Vturb(0 is a turbulent velocity at the scale I, k is 
the Boltzmann constant, fi = 0.6 is the mean molecular 



weight in a.m.u. 



is the proton mass, and T is the 



gas temperature. Both simulations and observational ev- 
idence suggest that, in regions around the virial radius, 
Pturb/^tiicrm ~ 0.1 ^ 0.3 at the t urbulence injection length 
scale lo of a few times 100 kpc fe.g.. lVazza et all201lf ). In the 
following, the value Pturb /ftherm = 0.2 will be used, which 
corresponds to a turbulent velocity i'turb(^o) = 640 km s~^ 
for a typical ICM temperature of 5 x 10^ K. 

For merger shocks, we wish to study scales Is smaller 
than the injection scale Iq. We set la/lo = 0.1, so that la is of 
the order of a few tens kpc, comparable to the relic width. 
The turbulent velocity «turb('s) is derived from Wturb(^o) by 
assuming Kolmogorov scaling, i.e. v (0 cx Z^/^ In this way. 



Slu 



«turb(0 _ Uturb(/o) 



1 D 



1 + M D 



1/3 



(6) 



The value of (5cj, for U/lo = 0.1 and i'turb('o) ~ 640 km s^^, 
is below 10~^^ s~^. For comparison, typical values of small- 
scale vorti city in turbulent regions of the IC M are about 
10"^'' s"^ (jKang et al.ll2007l : IPaul et al.ll201ll ;). For this rea- 
son, we can neglect the injection of vorticity due to shock 
curvature. Locally, the curvature term can be significant 
at locations where the derivative dMs/dS is large, namely 
where Ms varies a lot along the shock surface. An exam- 
ple is given by the large vorticity amplification at the inter- 
face between la rge-scale structu re filaments and propagating 
merger shocks (jPaul et al.i201ll ). Interestingly, bright notch- 



Figure 1. Comparison of the numerical coefBcients in front of the 
compressive and baroclinic terms (equation [TJ , described in the 
text, as a function of the shock Mach number Ms. The barocHnic 
coefficient has been multiplied by the turbulence Mach number 
Mt, whose value has been set to 0.3 (dashed line). The solid line 
refers to the compressive coefficient, namely the normalised den- 
sity ratio 11. 



like features have also been observed at the edges of some 
radio relics (|Rottgering et al.lll997l : iBagchi et al.ll2006l '). 

As for the baroclinic term, a qualitative estimate of the 
relative weight of this term as a function of the shock Mach 
number Ms, can be performed by a comparison of the nu- 
merical coefficients of the baroclinic and compressive terms 
in equation ([T]), which are 1/Ms (pi/(l-|-/i)M| — 1) and /i, re- 
spectively. Although a detailed analysis would require more 
information on the upstream flow and a complete compu- 
tation of the terms in equation ([T|, this task is beyond the 
scope of this paper. For consistency, we further multiply the 
numerical factor of the baroclinic term by Mt, in order to 
roughly take into account the magnitude of u in the product 
ujxu'm equation ([1} , where the flow velocity u is normalised 
to the sound speed. For the time being, 9(1/2 Mt)/dS in the 
baroclinic term is neglected. The normalised density ratio /i 
is expressed as a funct ion of Ms according to the Rankine- 
Hugoniot relations (cf. iKevlahan fc Pudrit j|2009l ): 



Mi 



1 + 1/2(7 - mi 



(7) 



The range of expected turbulent pressure ratios in clusters 
corresponds to a range of Mt between 0.2 and 0.35 on a 
scale of a few times 10 kpc. A value of 0.3 has been used 
for the comparison between the baroclinic and compressive 
coefficients of equation ([1} , shown in Fig. [T] 

As expected, p tends to 3 for increasing Ms- The am- 
plification of uj is mostly compressional. The baroclinic con- 
trib ution is small, as expected for a subsonic turbulent fiow 
(cf. iDel Sordo fc Brandenbur3l201l[ ). but not negligible, es- 
pecially at high values of Ms: for Mt = 0.3 and Ms = 5, 
the baroclinic coefficient is 40% of the compressional one. 
In this regime, both the compressional and baroclinic term 
contribute to the shock amplification of vorticity. 

In general, merger shocks are efficient in amplifying pre- 
existing vorticity, while they are less efficient in producing it 
via baroclinic effects. The estimate of the magnetic field in 
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the post-shock region thus depends on the upstream field, 
which we will estimate in the next Section. 



in the cluster periphery, a larger share of compressionally- 
driven turbulence might result in a smaller value for / in 
equation 



2.2 Upstream turbulent flow in the outer ICM 
and resulting magnetic field 

As known from MHD cosmological simulations, under rea- 
sonable assumptions for the seed field, it is not possible to 
reproduce the magnetic field expected in clusters only by 
adiabatic c ompression cause d by collapse during structure 
formation (|Dolag et al.ll20o"5h . In this Section, with the use 
of simple analytical estimates, we demonstrate that the tur- 
bulent dynamo is able to boost the level of ICM magnetisa- 
tion above the adiabatic compression values. These results 
are in rough agreement to the findings of MHD simulations 
(cf. Introduction). 

Turbulence is driven by structure formation at large 
length scales (of the order of Zq, a few times 100 kpc), and 
cascades to the dissipation length scale, presumably located 
at kpc or sub-kpc scales. All length scales in this range con- 
tribute to the amplification of B, but in a Kolmogorov cas- 
cade (or any other similar scaling) most of the kinetic energy 
is contained at the largest scal es, although the amp lification 
is faster for smaller eddies fcf. ISubramanian et al.ll2006 V In 
this simple analysis, we will restrict ourselves to the ampli- 
fication by the largest scale eddies {I — lo). 

The maximum magnetic field attainable by the tur- 
bulent dynamo is some fraction / of t he equiparti- 
tion v alue with the kin e tic en e rgy (see SubramanianI 



19981: ISubramanian et al] l2006l : ISchleicher et al.l l2010l : 



Federrath et~al]|201lD 



Bb. 



V/87rpekin = \//47rp«2 (8) 

In llapichino et al.l (|201ll l it is shown that the stirring in 
the cosmic gas is driven by a combination of compres- 
sional and solenoidal modes, with the latter ones dominat- 
ing in the ICM. For sub sonic, solenoidally driven turbulence 
iFederrath eTsU (|201ll ) derived / ~ 0.4 ; we will use this 
value in equation ((8|. 

Furthermore, we note that turbulence in clusters 

is expected not to be volume-filling. According to 

ISubramanian et al.l (|2006l l. for a volume-filling factor /v the 
rms magnetic field is equal to 



Bb, 



Bsat X /^^ 



V/v/47rpt;2 (9) 

Here we set /v = 0.3 (jiapichino fc Niemevejl2008l ). 

We use equation ([9]) for obtaining an estimate of the 
magnetic field in the ICM, before the shock propagation. 
We set Uturb(^o) ~ 640 km s~^, as in Section l2.ll As for 
the density in equation ((9|, the two values p = 10~^^ and 
lO"'^* g cm~'^ will be taken as representatives of the den- 
sity at central distances r = 0.5 iZvir and r = Rvir, respec- 
tively. The results are values of the saturation magnetic field 
Bsat.v = 2.5 fj,G at r = 0.5 Rvir, and 0.8 jj.G at r = 1 Rvir- 

These values are about a factor of two larger than the 
typical results from observations and MHD simulations, but 
our simple derivation cannot represent faithfully the whole 
complexity of the turbulent dynamo. Firstly, it is not guar- 
anteed that the field is saturated up to the scale Iq, and this 
could slightly reduce our estimates. Furthermore, especially 



2.3 Magnetic field and polarisation in the 
post-shock region 

Having shown that the compressional term is dominant for 
vorticity production across the shock, we go back to the 
question of the amplification of B, drawing some estimates 
from equations ^ and ((3]). Assuming that the upstream 
flow is isotropic, it is easy to see that the two B components 
parallel to the shock surface are amplified by a factor R, 
while the normal component is not amplified. The strength 
of the downstream field is thus 



2i?2 + 1 



(10) 



For R = 3, the shock amplification predicted by equa- 
tion p0|) is about 2.5, without taking into account some 
smaller contribution by baroclinicity. Given this amplifica- 
tion and the upstream estimates for B at the end of Section 
12.21 B reaches a post-shock value of 2.0 fj,G at r = -Rvir, 
and of 6.2 /iG at r = 0.5 -Rvir. For these estimates the same 
caveat mentioned at the end of Section [2. 21 applies. We defer 
a more detailed comment to Section [S] 

Another efi'ect of the compressional amplification is to 
make the downstream magnetic field anisotropic, by acting 
preferentially on the field component parallel to the shock. 
The degree of order introduced in the downstream field can 
be related to the polarisat i on ob served in radio relics. Ac- 
cording to iBonafede et all (|2009l ) (see also lBurnlll966l ). in 
this case the observed polarisation is 



3(5 3 



3S + 7 1 + {B?/Bi 



(11) 



where 5 — 2a + 1, a being the relic spectral index. The ran- 
dom and ordered components of the downstream magnetic 
fields are Bi and Bo, respectively. In the following, we label 
the components of B parallel to the shock surface as ordered, 
and the normal component as random; from equation (|10p . 
the ratio B^ / B^ is set therefore to 1/(2_R^). The first factor 
of the right-hand side of equation Ullf) can be interpreted as 
the intrinsic relic polarisation; for a typical value of a = 1, 
it is equal to 0.75. The second factor of the right-hand side 
is always smaller than unity and represents a suppression of 
the intrinsic polarisation. For our fiducial choice i? = 3, this 
factor is equal to 0.95, and Pobs = 0.71. 

In this oversimplified estimate we neglected to intro- 
ducc th e relic view ing angle, which further suppresses Pobs 
(|Enfilin et al.|[l99^ ). For the case of relics seen almost edge- 
on on the plane of sky, this suppression is not dominant, 
and in fact the polarisation derived from equation (lllf) is 
not far from that measured in CIZA J2242.8-f 5301 (50-60% 
level, [varLWceren ct al. 2010). In relics with a less favourable 
viewing geometry, or more complicated morphology, the po- 
larisation is lower (10 to 30%). 

Despite of the large uncertainties of the estimate in 
equation (|11|) . our aim here was to show that the ordering 
introduced by compressional amplification has the potential 



© 2012 RAS, MNRAS 000, [TH9] 



Amplification by shocks in clusters 5 



of producing a large polarisation of the radio emission, at 
odds with turbulent magnetic amplification, which acts in 
the opposite di rection of depolarising the signal, as observ ed 
in radio halos (jMurgia et al.ll2004l : iBonafede et al.ll2011bl '). 

Finally, we note that the possibility of a turbulent 
and magnetised pre-shock medium is not in contradiction 
with the regular morphology of many radio relics. Hydrody- 
namical simula tions of shock-turbulence interactions (e.g., 
iLee et al] 1 19931 ) show that the shock morphology is mainly 
unaltered, as long as Ms is sufficiently large in comparison 
to Mt. 



dissipation scale is ver y uncertain, but expe cted to be close 
to the kiloparsec scale (jSunvaev et al.|[2003h . 

The (unsurprising) result of this analysis is that the 
flow in the post-shock region, resulting from a laminar flow 
upstream, does not show fully developed turbulence, mainly 
because the timescale fovoi is too short. Clearly, if the up- 
stream flow itself is turbulent, mere compressional ampli- 
fication may suffice as shown in Section 12.11 As a side 
note, the compressional driving of turbulence, provided by 
a propagating shock, is rather inefficien t in amplifying B 
(Federrath ct al...2dTll : Ijones et al.ll201ll ). 



3 TURBULENT AMPLIFICATION IN THE 
POST-SHOCK REGION 

Here we show that, if the upstream flow is not turbulent, 
the stirring induced by the propagation of a single merger 
shock is not able to produce fully developed turbulence in 
the downstream region, and hence turbulent ampliflcation 
is not efficient. 

In hydrodynamical simulations of forced turbulence 
re.g.. lSchmidt et al.|[2009l : lFederrath et al.ll2009l ). turbulence 
is fully developed after a few large-scale eddy turnover times 
teddy, giveu by teddy = L/v, where L is the integral length 
(the largest length scale at which the flow is stirred) and v 
is the typical velocity at the scale L. 

During the propagation of the merger shock through the 
cluster outskirts, the downstream region is stirred. Let tcvoi 
be defined as the typical evolutionary timescale of the post- 
shock region. Turbulence is considered to be fully developed 
if the ratio 

N = ^ (12) 

teddy 

is larger than a few. Here, we interpret tcvoi as the advection 
timescale of the downstream gas over the width of the radio 
relic, Lroiic, or the region where a substantial amplification 
of the magnetic field is required. Thus we write 

_ -^rclic ^turb(-t/cddy) ^-^g-j 
-^'cddy ^down 

where Udown is the downstream velocity in the shock ref- 
erence frame. For shocks with Mach numbers between 2 
and 5 and shock velocities vs ^ 1000 — 4000 km s~^, 
Wdown = 1000 km s~^ is a reasonable value. Let us assume a 
maximal relic width of 200 kpc and a Kolmogorov scaling, 
i.e. v{l) oc l^^^ with a turbulent velocity of iiturb(200 kpc) = 
100 km s~^ which is motivated by simulation s of shock 
prop agation through a laminar upstream region (|Paul et al.l 
l201ll n. The turbulent velocity that we need to consider here 
is the turbulent velocity that results from a largely laminar 
infiow. Hence, we find that turbulence is not fully developed, 
i.e. A'^ < 2 as long as Leddy > 2 kpc. In cluster outskirts, the 

^ Tlie adaptive mesli refinement in lPauI et al ] l|201ll ') is designed 
to effectively refine the post-shock region, but it is not equally 
effective upstream. For this reason, h ere we cons i der th e down- 
stream turbulent velocity reported by I Paul et al. I ll201ll ') as rep- 
resentative of the stirring produced by a shock propagating in a 
nearly laminar upstream flow. This is exactly the case that wc 
want to study in this Section, although in real cluster environ- 
ments the turbulent velocity is expected to be much higher. 



4 DENSITY INHOMOGENEITIES IN THE 
UPSTREAM MEDIUM 

Based on Suzaku observations of the regions close to 
the virial radius, it has been suggested t hat there is 
some level of clumping in t he outer ICM jUrban et al.l 
imioncscu e t al.ri201lh . Numerical simulations by 
iNagai fc La u (20_11) have shown that clumping in the ICM, 
as for the substructure resolved in that work, starts being 
relevant at r ~ R200, and thus affects mostly the outer 
part of the distance range studied throughout this paper. 
Gas clumps in th e ICM have recently been invoked by 
iBirnb 01m fc Dekell ([ 201 ih in their mechanism for balancing 
radiative cooling in cool-core clusters. Moreover, they esti- 
mated that about 5 per cent of the accreted baryons are in 
gas clumps with masses of ~ 10* Mq. 

As well as for the turbulent medium, the propagation of 
a shock through an inhomogeneous medium can be a source 
of additional vorticity and thus magnetic amplification. Vor- 
ticity is generated by the baroclinic mechanism at de nsity in- 
homogeneities (Samtancy & Zabusky 1994). In Ino ue et al.l 
(l2009h . for conditions typical of supernova remnants, the 
interaction between a strong shock wave and a multiphase 
interstellar medium produces amplifications well above the 
level attainable by turbulence. 

Here, we make a few geometrical estimates: a sphere 
with R = 1 kpc and density pb ~ XP, where the density 
contrast x is set to 10 and p is the ambient density in the 
outer ICM (set to 10"^* g cm"^, i.e. about 200 times the 
average baryon density), has a mass Mb = 6 X 10* Mq. As 
a working hypothesis, we assume that 5% of the baryons 
in the outer ICM is in such clouds. We consider a vol- 
ume of (0.05 X 1.7 X 1.7) Mpc, similar to the volume of 
the CIZA J2242. 8-1-5301 relic, with an average density p. 
Some simple algebra shows that in this volume there are 
1.8 X 10^ such clumps. From geometrical considerations, 
their average spac i ng wi ll be around 10 kpc. According to 
iPoludnenko et ahl l|2002l '). the ratio between the clump ra- 
dius and this typical separation is smaller than the minimum 
value needed for the clumps to interact and merge in their 
late evolution. The initial volume fraction occupied by the 
clumps is Vb.in = 0.005. By visual inspection of the results of 
IPoludnenko et al.l l|2002l l. it is assumed that the clump gets 
spread in a volume Vb,flnai ~ lOVb.in before being dispersed. 
This results in a final volume filling factor which remains 
small (below the 10% level). 

By fiux conservation, we set the initial field in the 
clumps as -Bb,in = Bup{pb/p)'^^^ ~ 3-Bup, where Bup is the 
upstream magnetic field in the ICM. 



© 2012 RAS, MNRAS 000, [TH9] 



6 L. lapichino and M. Briiggen 



In the i r MH D simulations of shock-cloud interaction, 
ih in et al.l (|2008t ) find that, after the shock propagation and 
before the complete clump dispersal, the clump magnetic 
field is amplified up to a value Bb.Hnai ~ /bSb.in, with the 
amplification factor /b which depends on the orientation of 
B and on the /3 parameter, defined as /3 = P/Pb, the ratio 
between the thermodynamical and magnetic pressure. For 
the problem parameters considered here. 



P _ Stt pb fcTb 
Pb tJ-mpBl^^ 



(14) 



where we used Tb = 5 x lO'^ K (smaller than the ambient T 
by a factor of x for ensuring pressure equilibrium between 
the ICM and the clumps), and -Bb.in ~ 1-5 fJ,G, resulting 
from Bup = 0.5 /iG. The value of (5 derived i n equ ation 
l|14p is within the range examined bv lShin et al.l (|2008l ): the 
corresponding amplification factor is set to /b = 2 (in line 
with the cited res ults, but somewh at extreme). 

According to lShin et all ||2008D and many other similar 
studies, the typical timescale for the clump disruption (and 
for the B amplification) is set by the so-called cloud-crushing 
time, which is equal to clump shock crossing timescale times 
X^^^- Given the large difi'erence in size between the clump 
and the post-shock region, it is obvious that the dispersal 
mechanism (and the related B amplification) has enough 
time to develop during the shock propagation in the post- 
shock region. 

Summarising the arguments provided above, we arrange 
them in an estimate of the B amplification. First, we set 
Bdown = 2.5-Bup as the downstream value of amplified the 
magnetic field, without the effect of clumping (cf. equation 
nop . The volume- weighted downstream B, including clump- 
ing, is (1 — Vb,flnal)Bdown + Vb.flnal (-Bb, final ) ~ 1.07 -Bdown • 

The whole effect of the clumping considered here is an am- 
plification at the negligible level of some percent. 

There are many arbitrary parameters in this analysis, 
but the amplification level will not be significantly larger 
without invoking a drastically higher baryon fraction in 
these clumps, or larger sizes. The latter is unlikely as it 
would destroy the smoo th morphology of the radio relic 
l|van Weeren et al.|[2oTlah . 



5 DISCUSSION AND CONCLUSIONS 

In this work, we examined the amplification of magnetic 
fields in the post-shock region of merger shock waves in 
galaxy clusters. We considered several mechanisms for am- 
plifying magnetic fields in the downstream region of the 
shock, including small-scale dynamo action, compressional 
amplification, baroclinic amplification and amplification by 
the interaction with density inhomogeneities. 

We demonstrated that turbulent amplification is ineffi- 
cient in producing magnetic fields of the magnitude required 
by observations in radio relics because the stirring induced 
by the shock propagation in the post-shock region is not 
sufficient to drive fully developed turbulence. Instead, we 
propose that the shock propagates through regions of the 
outer ICM with a turbulent pressure support of the order of 
10 to 30 percent. 

The dominant mechanism is the compressional amplifi- 
cation, which amplifies the magnetic field components par- 



allel to the shock surface. Together with compression at 
shocks, baroclinic effects are expected to be marginally rel- 
evant for high Mach numbers (Ms > 3) and intense ICM 
turbulence (Mt = 0.3). Motivated by recent observational 
indications, we explored in Section |4] the role of clumping 
for the magnetic field amplification, which turns out to be 
modest. 

Compressional amplification can suffice to explain mag- 
netic fields in radio relics. Upstream fields are amplified by 
the outgoing merger shock, as expressed in Section [2] by 
equation (|10|l . For a typical shock Mach number AIs = 3, 
the resulting amplification factor is equal to 2.5. The mag- 
netic field in the post-shock region is thus estimated to be 
6.2 ^G at r = 0.5 i?vir, and 2.0 /iG at r = 1 iivir. 

These values are in good agreement with the avail- 
able estimates of magnetic fields in radio relics (e.g., 
iBonafede et all l2009l : Ivan Weeren et aLlbOllbl ). In the case 
of CIZA J2242. 8-1-5301, the field from observations is be- 
tween 5 and 7 fiG, which is well matched by our predic- 
tion. There is some discrepancy between observations and 
theory for the NW relic of A3667, where the observational 
estimate for the relic fi eld is B = 3 fj,G at r ~ 1 7?vir 
(|Finoguenov et al]l2010l ). to be compared to a theoretical 
prediction at the level of about 2 ^G. Given the simplicity 
of our model, this still constitutes reasonable agreement. 

The proposed model for the compressional amplifica- 
tion of the magnetic field, in contrast with turbulent dy- 
namo models, can naturally account for a large polarisation 
of the radio emission (equation llip , although more realis- 
tic est imates should incor porate the role of the relic viewing 
angle l|Enfilin et aLlliggsl ). 

Compared to X-ray observations, the outer regions 
of clusters can potentially be better studied by Sunyaev- 
Zel'dovich (SZ) observations because the SZ effect depends 
only on the electron density to the first power, while the X- 
ray emission depends on the electron density squared. Fa- 
cilities such as the Atacama Cosmology Telescope (ACT), 
the South Pole Telescope (SPT) and the Planck satellite are 
searching for the SZ signal of galaxy clu sters and have foun d 
some interesting first results (see, e.g.. iBagchi et aLll201ll '). 
For example, it was found that the SZ signal caused by clus- 
ters is by a factor of 2 smaller than predicted by mod- 
els of clusters that suggest th at the pressure by the elec- 
trons has been overpredicted l|Lueker et al.|[201oh . This is 
presumably caused by a substantial nonthermal pressure at 
cluster outski rts, most of it is likely to be turbulent pres- 
sure (see e.g. IShaw et"al] |2010(). This picture is also sup- 
ported by simul ations ( Lau et akl 120091 : IVazza et al. I I2OO9I : 
iPaul et al]|201ll ). However, the volume-filling fa ctor of tur- 
bulence in the cluster outs kirts is still debated l|Valdarninil 
I2OIII : llapichino et"al]|201ll ) . 

Currently, there are few measurements of magnetic 
fields in regions so far from the cluster centre, such as for 
example I Clarke et al.l l|200ll ). who detected some Faraday 
RM excess out to central distances around 0.5 Mpc . 
However, with the advent of the SKA, the magnetic fields 
in such regions can be better probed, by measuring the 
Faraday rotation in a much larger number of background 
sources (|Krause et al.l [20091 ). Synchrotron radiation propa- 
gating through a magnetised plasma undergoes Faraday ro- 
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tation, and the rotation measure RM is given by: 



RM = 812 



ncBudl rad m 



(15) 



where Lpath is the path along the line of sight in units of kpc, 
ric = p/{mpfic) is the electron density in cm~^ {rrip is the 
proton mass, — 1-14 is the average mass per electron in 
a.m.u), and -By is the component of B along the line of sight, 
in fiG. We assume that a tangled magnetic field topology 
in the ICM leads to a Gaussian distribution of RMs with 
a vanishing mean ( RM) and a RM variance <t^m given by 
l|Murgia et al.ll2004i ): 



2 

CRM 



(RM^) 



(ne-B||)^d/ 



(16) 



where A is the typical size of a turbulent cell. Under the 
assumption that the gas density profile follows a /3-model 



/, , 2 , 2x3/3/2 

ric = no (l + r /r^) 



(17) 



with central density no an d core radius rc , the R M dispersion 
ctrm can be expressed as (|Fehenlll99d see also iDolag et alj 
I2OOII ): 



0"RM 



KBnory^A^^^ / r(3/3-0.5) 

(i+^2/^2)(6^-i)/4Y rm 



(18) 



where F is the Gamma function, and K is a, numerical 
factor which depends on the location of the radio source 
with respect to the cluster, and is equal to 624 or 441 for 
background or embedded sources, respectively. The mag- 
netic field B is assumed to be constant in the ICM, al- 
though a radial depen dence of the type B{r) = Bo[n{r)/no]^ 
l|Bonafede et al.ll2010l ) can be i ntroduced by subs tituting /? 
with /3(1 + rj) in equation HH)) (iDolag et al1l200ll ). 

We adopt the sam e /?-model parameters as in 
Ivan Weeren et al.l (|2011al ): Tc = 134 kpc, /3 = 2/3, m = 
2 X 10~^ cm~^, and we set r = 1.5 Mpc. The typical size 
of the turbulent eddies A is set to 100 kpc, according to 
ISubramanian et al.l (|2006t ) and to the assumptions of this 
work. Finally, we consider a background source {K = 624). 
The resulting RM dispersion for the assumed upstream field 
of -B = 2.5 fiG is CTRM = 9, and for a shock-amplified field 
of B = 6.2 fiG, it is ctrm = 22. 

This estimate is probably too rough because the mag- 
netic field has the amplified value only in the post-shock 
region, and not everywhere in the ICM along the line of 
sight. A further reduction of ctrm comes from imposing a 
radial profile for B: for a line of sight at r — 0.5 i?vir from 
the cluster centre, it results in ctrm about 20% smaller than 
in equation (|18[> . as es timated from applyin g it on Coma, 
using parameters from ISonafede et al.l l20ld . On the other 
hand, for a line of sight passing through a radio relic, the 
assumption of a Gaussian distribution for the RM might not 
hold, because of the ordering introduced on B by the com- 
pressional amplification. As a result, (RM) could be non- 
vanishing for lines of sight crossing the relic. 

The problem of Faraday rotation in connection to clus- 
ter merger shocks is interesting and worth being further 
studied, also by means of numerical simulations, because 
the derived values of ctrm, despite of uncertainties in their 
derivation, are within the sensitivit y goal of upcoming in- 
struments such as SKA or ASKAP jseck fc Gaensleill2003 : 
iKrause et alll2009l '). 



Besides the mechanisms considered here, there are a 
number of other processes that are able to amplify the mag- 
netic field at shock s. Viable channels are the Weibel insta- 
bility (|Weibellll959h and the amplification dri ven by cosmic- 
ray a cceleration (the so-called Bell instability; iLucek fc Belj 
I2OOOI : Bell 2004). The main issue with these mechanisms is, 
as in the case of clumping, the length scale: this is a known 
concern for the Weibel instability (Medv edev et al. 2006), 
which has the fastest growing mode below the parsec scale 
by orders of magnitude, although it has been speculated 
about a rapid transfer of ma gnetic energy to cosmological 
scales l|Medvedev et al.]|2005l '). As for the Bell instability, in 
the non-linear stage the do minant length scale can grow u p 
to the CR Larmor radius (jRiguelme fc Spitkovskvl l2009l l. 
For the typical conditions considered in this work, we argue 
that this length scale r emains below or around the p arsec 
scale (cf. equation 6 of iRiguelrne fc Spitkovsk^l2009l . and 
the discussi on therein). 

Finally, iBeresnvak et al.l (|2009l ) studied small-scale dy- 
namos in the context of diffusive shock acceleration. They 
consider the effect of the CR shock precursor on the up- 
stream density perturbations, and identify a mechanism that 
is able to generate an efficient amplification of the mag- 
netic field in the precursor. In this model, an important 
condition is the existence of density perturbations in the 
upstream flow. In the framework of ISM studies, this is 
guaranteed by the relatively large Mach number of the fiow. 
However, in subsonic flow, as in the ICM, the amplitude 
of the density perturbation s is much smaller (see for ex- 
ample Fedcrrath ct al. 2008), and the amplification in the 
precursor inefficient. For length scales and separations, the 
clumping described in Sect ion [4] cannot contribut e to this 
mechanism. As a remedy, IBeresnvak et al.l (|2009l ') hint at 
the possibility that density fluctuations are produced in the 
CR precursor itself. This needs to be studied in future work. 

Faraday rotation measurements in the outer ICM will 
help to discriminate between purely hydrodynamical models 
of field amplification, which require some level of turbulence 
and magnetisation in the cluster outskirts, and other models, 
which link magnetic field amplification to CR acceleration. 



ACKNOWLED CEMENTS 

MB acknowledges support from the Deutsche Forschungs- 
gemeinschaft from the grant FOR1254. We are grateful to 
the referee, Klaus Dolag, for his insightful and construc- 
tive suggestions, which improved the presentation of this 
work. Thanks to V. Antonuccio-Delogu and J. C. Niemeyer 
for having supported the first phases of this work, to 
A. Bonanno and W. Schmidt for useful comments on the 
manuscript, and to M. Bartelmann, C. Federrath, R. Klessen 
and F. Vazza for interesting discussions. 



REFERENCES 

Bagchi J., Durret F., Neto G. B. L., Paul S., 2006, Science, 
314, 791 

Bagchi J., Sirothia S. K., Werner N., Pandge M. B., Kan- 
tharia N. G., Ishwara-Chandra C. H., Gopal-Krishna, 
Paul S., Joshi S., 2011, ApJ, 736, L8 



© 2012 RAS, MNRAS 000, [TH9] 



8 L. lapichino and M. Briiggen 



Battaglia N., Pfrommer C, Sievers J. L., Bond J. R., 

Enfilin T. A., 2009, MNRAS, 393, 1073 
Beck R., Gaensler B. M., 2004, New Astronomy Review, 

48, 1289 

Bell A. R., 2004, MNRAS, 353, 550 

Beresnyak A., Jones T. W., Lazarian A., 2009, ApJ, 707, 
1541 

Birnboim Y., Dekel A., 2011, MNRAS, 415, 2566 
Bonafede A., Dolag K., Stasyszyn F., Murante G., Borgani 

S., 2011a, MNRAS, 418, 2234 
Bonafede A., Feretti L., Murgia M., Govoni F., Giovannini 

G., Dallacasa D., Dolag K., Taylor G. B., 2010, A&A, 513, 

A30 

Bonafede A., Giovannini G., Feretti L., Govoni F., Murgia 

M., 2009, A&A, 494, 429 
Bonafede A., Govoni F., Feretti L., Murgia M., Giovannini 

G., Briiggen M., 2011b, A&A, 530, A24 
Briiggen M., Ruszkowski M., Simionescu A., Hoeft M., 

Dalla Vecchia C., 2005, ApJ, 631, L21 
Burn B. J., 1966, MNRAS, 133, 67 

Clarke T. E., Kronberg R P., Bohringer H., 2001, ApJ, 547, 
LIU 

Davies G., Widrow L. M., 2000, ApJ, 540, 755 
Del Sordo F., Brandenburg A., 2011, A&A, 528, A145 
Dolag K., Bartelmann M., Lesch H., 2002, A&A, 387, 383 
Dolag K., Grasso D., Springel V., Tkachev I., 2005, Journal 

of Cosmology and Astro-Particle Physics, 1, 9 
Dolag K., Schindler S., Govoni F., Feretti L., 2001, A&A, 

378, 777 

Dolag K., Stasyszyn F., 2009, MNRAS, 398, 1678 
Donnert J., Dolag K., Lesch H., Miiller E., 2009, MNRAS, 
392, 1008 

Dubois Y., Teyssier R., 2008, A&A, 482, L13 

Enfilin T. A., Biermann P. L., Klein U., Kohle S., 1998, 

A&A, 332, 395 
Enfilin T. A., Briiggen M., 2002, MNRAS, 331, 1011 
Enfilin T. A., Gopal-Krishna, 2001, A&A, 366, 26 
Federrath C, Chabrier G., Schober J., Banerjee R., Klessen 

R. S., Schleicher D. R. G., 2011, Phys. Rev. Lett., 107, 

114504 

Federrath C., Klessen R. S., Schmidt W., 2008, ApJ, 688, 
L79 

— , 2009, ApJ, 692, 364 

Felten J. E., 1996, in Astronomical Society of the Pacific 
Conference Series, Vol. 88, Clusters, Lensing, and the Fu- 
ture of the Universe, V. Trimble & A. Reisenegger, ed., p. 
271 

Ferrari C, Govoni F., Schindler S., Bykov A. M., Rephaeli 

Y., 2008, Space Sci. Rev., 134, 93 
Finoguenov A., Sarazin C. L., Nakazawa K., Wik D. R., 

Clarke T. E., 2010, ApJ, 715, 1143 
Giacintucci S., Venturi T., Macario G., Dallacasa D., 

Brunetti G., Markevitch M., Cassano R., Bardelli S., 

Athreya R., 2008, A&A, 486, 347 
Govoni F., Murgia M., Feretti L., Giovannini G., Dolag K., 

Taylor G. B., 2006, A&A, 460, 425 
Hoeft M., Briiggen M., 2007, MNRAS, 375, 77 
Hoeft M., Briiggen M., Yepes G., Gottlober S., Schwope 

A., 2008, MNRAS, 391, 1511 
lapichino L., Niemeyer J. C, 2008, MNRAS, 388, 1089 
lapichino L., Schmidt W., Niemeyer J. C, Merklein J., 

2011, MNRAS, 414, 2297 



Inoue T., Yamazaki R., Inutsuka S., 2009, ApJ, 695, 825 
Jones T. W., Porter D. H., Ryu D., Cho J., 2011, ArXiv 

e-prints, 1108.1369 
Kang H., Ryu D., 2011, ApJ, 734, 18 

Kang H., Ryu D., Cen R., Ostriker J. P., 2007, ApJ, 669, 
729 

Kempner J. C, Blanton E. L., Clarke T. E., Enfihn T. A., 
Johnston-HoUitt M., Rudnick L., 2004, in The Riddle 
of Cooling Flows in Galaxies and Clusters of galaxies, 
T. Reiprich, J. Kempner, & N. Soker, ed., p. 335 
Kevlahan N., 1997, Journal of Fluid Mechanics, 341, 371 
Kevlahan N., Pudritz R. E., 2009, ApJ, 702, 39 
Krause M., Alexander P., Bolton R., Geisbiisch J., Green 

D. A., Riley J., 2009, MNRAS, 400, 646 
Lau E. T., Kravtsov A. V., Nagai D., 2009, ApJ, 705, 1129 
Lee S., Lele S. K., Moin P., 1993, Journal of Fluid Mechan- 
ics, 251, 533 
Lucek S. G., BeU A. R., 2000, MNRAS, 314, 65 
Lueker M., Reichardt C. L., Schaffer K. K., Zahn O., Ade 
P. A. R., Aird K. A., Benson B. A., Bleem L. E., Carlstrom 
J. E., Chang C. L., Cho H.-M., Crawford T. M., Crites 
A. T., de Haan T., Dobbs M. A., George E. M., Hall N. R., 
Halverson N. W., Holder G. P., Holzapfel W. L., Hrubes 
J. D., Joy M., Keisler R., Knox L., Lee A. T., Leitch E. M., 
McMahon J. J., Mehl J., Meyer S. S., Mohr J. J., Montroy 
T. E., Padin S., Plagge T., Pryke C, Ruhl J. E., Shaw L., 
Shirokoff E., Spieler H. G., Stalder B., Staniszewski Z., 
Stark A. A., Vanderlinde K., Vieira J. D., Williamson R., 

2010, ApJ, 719, 1045 

Macario G., Markevitch M., Giacintucci S., Brunetti G., 
Venturi T., Murray S. S., 2011, ApJ, 728, 82 

Markevitch M., Govoni F., Brunetti G., Jerius D., 2005, 
ApJ, 627, 733 

Medvedev M. V., Fiore M., Fonseca R. A., Silva L. O., Mori 

W. B., 2005, ApJ, 618, L75 
Medvedev M. V., Silva L. O., Kamionkowski M., 2006, ApJ, 

642, LI 

Murgia M., Govoni F., Feretti L., Giovannini G., Dallacasa 
D., Fanti R., Taylor G. B., Dolag K., 2004, A&A, 424, 429 

Nagai D., Lau E. T., 2011, ApJ, 731, LIO 

Nuza S. E., Hoeft M., van Weeren R. J., Gottlober S., Yepes 
G., 2012, MNRAS, 420, 2006 

Paul S., lapichino L., Miniati F., Bagchi J., Mannheim K., 

2011, ApJ, 726, 17 

Poludnenko A. Y., Frank A., Blackman E. G., 2002, ApJ, 
576, 832 

Riquelme M. A., Spitkovsky A., 2009, ApJ, 694, 626 
Rottgering H. J. A., Wieringa M. H., Hunstead R. W., 

Ekers R. D., 1997, MNRAS, 290, 577 
Ryu D., Kang H., Cho J., Das S., 2008, Science, 320, 909 
Samtaney R., Zabusky N. J., 1994, Journal of Fluid Me- 
chanics, 269, 45 
Schleicher D. R. G., Banerjee R., Sur S., Arshakian T. G., 
Klessen R. S., Beck R., Spaans M., 2010, A&A, 522, A115 
Schmidt W., Federrath C, Hupp M., Kern S., Niemeyer 

J. C, 2009, A&A, 494, 127 
Shaw L. D., Nagai D., Bhattacharya S., Lau E. T., 2010, 

ApJ, 725, 1452 
Shin M.-S., Stone J. M., Snyder G. F., 2008, ApJ, 680, 336 
Simionescu A., Allen S. W., Mantz A., Werner N., Takei Y., 
Morris R. G., Fabian A. C, Sanders J. S., Nulsen P. E. J., 
George M. R., Taylor G. B., 2011, Science, 331, 1576 



© 2012 RAS, MNRAS 000, [TH9] 



SkiUman S. W., Hallman E. J., O'Shea B. W., Bums J. O., 
Smith B. D., Turk M. J., 2011, ApJ, 735, 96 

Subramanian K., 1998, MNRAS, 294, 718 

Subramanian K., Shukurov A., Haugen N. E. L., 2006, MN- 
RAS, 366, 1437 

Sunyaev R. A., Norman M. L., Bryan G. L., 2003, Astron- 
omy Letters, 29, 783 

Truesdell C, 1952, J. Aero. Sci., 19, 826 

Urban O., Werner N., Simionescu A., Allen S. W., 
Bohringer H., 2011, MNRAS, 414, 2101 

Valdarnini R., 2011, A&A, 526, A158 

van Weeren R. J., Briiggen M., Rottgering H. J. A., Hoeft 

M., 2011a, MNRAS, 418, 230 
van Weeren R. J., Hoeft M., Rottgering H. J. A., Briiggen 

M., Intema H. T., van Velzen S., 2011b, A&A, 528, A38 
van Weeren R. J., Rottgering H. J. A., Bagchi J., Ray- 

chaudhury S., Intema H. T., Miniati F., Enfilin T. A., 

Markevitch M., Erben T., 2009, A&A, 506, 1083 
van Weeren R. J., Rottgering H. J. A., Briiggen M., Hoeft 

M., 2010, Science, 330, 347 
Vazza F., Briiggen M., van Weeren R., Bonafede A., Dolag 

K., Brunetti G., 2012, MNRAS, 2652 
Vazza F., Brunetti G., Gheller C., Brunino R., Briiggen 

M., 2011, A&A, 529, A17 
Vazza F., Brunetti G., Kritsuk A., Wagner R., Gheller C., 

Norman M., 2009, A&A, 504, 33 
Weibel E. S., 1959, Physical Review Letters, 2, 83 
Xu H., Li H., Collins D. C., Li S., Norman M. L., 2010, 

ApJ, 725, 2152 
— , 2011, ApJ, 739, 77 

This paper has been typeset from a T^/ I^TJ^X file prepared 
by the author. 



Amplification by shocks in clusters 9 



© 2012 RAS, MNRAS 000, [THH] 



